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i Ma p EE Hb Voronoi 
子 块 体 单元 DDA 方法 模拟 
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摘 要 : 非 连续 变形 分 析 (DDA ) 方法 是 计算 离散 可 变形 块 体系 统 力 学 行为 的 数值 计算 方法 ,可 通过 
子 块 体 单 元 DDA 方法 模拟 岩石 的 开裂 破坏 。 考 虑 到 Voronoi 多 边 形 颗粒 与 细 观 尺度 下 岩石 矿物 蝇 
粒 形 态 的 相似 性 ,提出 一 种 基于 随机 圆 的 Voronoi 颗粒 单元 模型 生成 方法 ;并 通过 完整 及 带 预 制 列 
纹 岩 石 圆 盘 径 向 压缩 破坏 的 模拟 ,验证 岩石 破裂 问题 Voronoi 子 块 体 单元 DDA 模拟 方法 的 适用 性 。 
结果 表明 , 当 子 块 体 单元 数 较 小 时 , 圆 盘 表现 出 更 高 的 整体 强度 ; 随 着 子 块 体 单元 数 的 增 大 ,起 裂 处 
位 置 更 接近 真实 ,开裂 破坏 路 径 更 清晰 ; 子 块 体 单元 数 较 大 时 不 同 倾角 预制 裂纹 圆 盘 破坏 的 模拟 结 
* 果 与 实验 结果 高 度 吻合 ,并 能 有 效 反 映 圆 盘 中 心 加 工 小 孔 对 开裂 破坏 路 径 的 影响 。 使 用 Voronoi F 
块 体 单元 DDA 方法 能 够 有 效 模 拟 岩 石 的 开裂 破坏 过 程 , 为 进一步 开展 基于 Voronoi 颗粒 单元 模型 
的 岩石 开裂 破坏 模拟 创造 了 条 件 。 

关键 词 :岩石 圆 盘 ;开裂 破坏 ;Voronoi 颗粒 单元 ; 非 连 续 变形 分 析 (DDA ) 
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Simulation of rock disc failure under diametrical compression 
by the Voronoi sub-block element DDA 


ZHOU Dabo! , NING Youjun’ ,LI Chunyu , WANG Daguo? 


(1. School of Mechatronic Engineering , Southwest Petroleum University ,610500 Chengdu , China ; 


2. School of Civil Engineering and Geomatics , Southwest Petroleum University ,610500 Chengdu , China) 


Abstract ; Discontinuous deformation analysis ( DDA ) is a numerical method that computes the mechanical 
behaviors of discrete deformable-block systems , which can be used to simulate the cracking failure of rock 
by the sub-block method. Considering the similarity of morphology between Voronoi polygonal particles and 
rock mineral grains shapes at the mesoscale ,a generation method of Voronoi particle element model based 
on random circles is proposed. By simulating the cracking failure process of intact rock disc and rock discs 
with a pre-existing crack , the feasibility of simulating rock cracking failure with the Voronoi sub-block ele- 
ment DDA is verified. Results indicate that the rock disc shows higher integral strength with the decrease 


of the number of sub-block elements. With the increase of the element number ,the crack initiation position 
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is closer to that in reality , and the failure path gets clearer. With a large-enough number of the elements, 


the cracking failure simulation results of the discs with different pre-existing crack inclinations are highly 


consistent with the experimental results , and the influence of the center machining hole in the disc on the 


failure path can be accurately captured as well. The Voronoi sub-block element DDA can effectively simu- 


late the cracking failure process of rock , which provides conditions for the further simulation of rock crack- 


ing failure with Voronoi particle elements. 


Key words rock discs ;cracking failure; Voronoi particle element ;discontinuous deformation analysis (DDA) 


岩石 是 一 种 不 连续 性 、 非 均 质 、 各 向 异性 的 天 然 
材料 ,岩石 在 外 载 作用 下 的 破坏 过 程 实质 上 是 岩石 
中 裂纹 的 萌生 、 发 育 与 贯通 的 过 程 。 圆 盘 劈 裂 是 典 
型 的 岩石 变形 破 袭 实验 方法 5031 ,而 数值 模拟 是 研究 
岩石 变形 破裂 行为 的 重要 手段 。 

在 连续 数值 计算 方法 中 ,有 限 元 法 (finite ele- 
mánt method , FEM ) 在 岩石 破坏 问题 的 研究 中 得 到 了 
Vae; Hi. ZHU 等 “采用 基于 FEM 的 RFPA 方 
法 种 拟 了 岩石 圆 盘 在 静态 和 动态 加 载 下 的 破坏 过 
FÉSSAKSALA 等 "对 岩石 圆 盘 进 行 不 同 冲 击 速 度 的 
埃 愉 实验 并 开展 了 FEM 模拟 。FEM 中 基于 单元 损 
价 尖 变形 或 单元 删除 的 破坏 模拟 方法 不 是 对 岩石 破 
裂 面 的 直接 描述 ,对 模拟 后 续 的 岩石 大 位 移 及 接触 
问题 存在 困难 。 在 非 连 续 数 值 计算 方法 中 ,离散 元 
Vet discrete element method , DEM ) 在 和 岩石 的 开裂 破 
坏 模 拟 中 也 得 到 了 广泛 应 用 。 如 孟 京 京 等 使 用 
DEM 颗粒 流程 序 ( particle flow code, PFC) 模拟 了 不 
同 忠 心 角 平台 圆 盘 的 巴西 臂 裂 和 直接 拉 伸 试 验 ; 
YANG 等 中 利用 PFC 模拟 了 完整 及 带 预 制 裂纹 岩石 
圆 各 的 破裂 过 程 。PFC 将 岩石 材料 离散 成 圆 形 果 
粒 6 基 于 粒 间 的 接触 描述 材料 的 拉 、 前 破坏 ,存在 颗 
粒 形态 与 岩石 晶体 形态 差异 较 大 ,以 及 颗粒 尺寸 和 
接触 参数 较 难 确定 等 难题 。 

连续 - 非 连续 耦合 或 统一 的 数值 方法 对 于 岩石 
开裂 破坏 的 模拟 具有 较 好 的 适用 性 。 如 张 世 瑞 
等 中 基于 FEM-DEM 耦合 方法 研究 了 不 同 Voronoi 
颗粒 尺寸 对 岩石 试 样 单 轴 抗 压强 度 及 损伤 应 力 的 影 
响 。 流 形 元 法 (numerical manifold method ,NMM ) 是 
连续 和 非 连 续 相 统一 的 数值 方法 ,被 应 用 于 岩石 等 
材料 的 复杂 开裂 破坏 问题 模拟 522 。LIU Age 
用 NMM 法 模拟 了 岩石 的 单 轴 压 缩 与 圆 盘 辟 裂 等 经 
HIS ,结果 表明 采用 Voronoi 颗粒 能 够 较 好 地 反映 
矿物 粒度 对 岩石 力学 性 能 的 影响 。 

韭 连 续 变形 分 析 ( discontinuous deformation anal- 


ysis, DDA) ”是 能 够 模拟 任意 形状 块 体 单 元 的 隐 式 


=H 


非 连续 数值 计算 方法 , 较 多 学 者 采用 DDA 方法 模拟 
Tun m dE sk e UU V n JIAO 等 ' 通过 引 
和 人 虚拟 节理 的 概念 ,采用 三 角形 子 块 体 单元 ,模拟 了 
AUS Ré vt UU B) ESL Jp REDE. fü 
克 松 等 ”在 DDA 子 块 体 单元 开裂 算法 中 根据 邻近 
子 块 体 的 应 力 状态 进行 开裂 判断 ,降低 了 网 格 依赖 
KE, YU 45! fg DDA 中 提出 了 一 种 分 布 式 弹簧 接 
触 模型 ,提升 了 接触 强度 问题 的 模拟 精度 。 张 开 雨 
等 ”采用 基于 Voronoi 离散 的 DDA 方法 研究 了 脆 
性 材料 的 动态 裂纹 扩展 。 

本 研究 发 展 一 种 基于 随机 圆 的 Voronoi JUR ES 
散 模型 建立 方法 ,并 开展 完整 和 带 预 制 裂纹 岩石 圆 
盘 径 向 压缩 破坏 的 DDA 模拟 ,研究 Voronoi 子 块 体 
单元 尺寸 对 模拟 结果 的 影响 ,验证 Voronoi 子 块 体 单 
元 DDA 方法 模拟 岩石 开裂 破坏 的 可 行 性 。 


1 Voronoi 颗粒 单元 模型 的 建立 方法 


图 1 为 Voronoi 颗粒 单元 模型 的 建立 过 程 。 步 
又 一 通过 在 子 域 中 生成 给 定 级 配 的 随机 圆 引 入 级 配 
特性 控制 参数 ,并 通过 子 域 的 复制 以 及 步 又 三 中 的 
中 心 区 域 截取 获得 模型 周期 特性 。Voronoi 颗粒 在 
步骤 二 中 基于 随机 圆 的 圆心 生成 ,进而 通过 步骤 三 
中 的 中 心 截取 区 域 再 拼接 获得 足够 大 的 Voronoi 颗 
粒 单元 分 布 区 域 。 


步 又 一 : 按 给 定 级 配 在 子 域 中 生成 随机 
列 并 进行 子 域 复 制 


随机 圆 圆心 生成 Voronoi 
JU 


步骤 三 : 截取 出 周期 性 区 域 并 进行 拼接 


: 根据 实际 问题 切割 生成 


E Voronoi 颗 粒 单元 模型 


步骤 五 :转化 生成 DDA 计 算 模 型 


图 1 Voronoi 颗粒 单元 模型 的 建立 过 程 


Fig.l Establishment process of Voronoi particle element model 
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与 直接 生成 全 域 的 Voronoi 颗粒 单元 相 比 ,拼接 
的 方法 具有 更 高 的 生成 效率 。 进 而 在 步骤 四 中 , 根 
据 实际 问题 的 需要 ,通过 切 制 生成 任意 形状 ( 含 带 孔 
洞 的 情况 ) 的 岩石 Voronoi 颗粒 单元 模型 。 最 后 ,在 
步骤 五 中 根据 DDA 块 体 单元 的 格式 ,将 Voronoi 颗 
粒 单 元 模型 转化 DDA 计算 模型 。 

以 图 2 为 例 ,在 步骤 一 中 ,首先 在 Al 正方 形 子 
域内 按 给 定 级 配 信 息 生成 互 不 重 和 到 的 随机 圆 ,并 将 
Al 子 域内 的 随机 圆 复制 到 相 邻 的 A2、A3 和 A4 F 
域 。 在 Al 子 域 随 机 圆 投放 的 过 程 中 ,实际 的 投放 区 
域 将 向 外 延伸 最 大 尺寸 圆 形 颗粒 直径 的 0.25 倍 , 以 
尽量 消除 子 域 边 界 处 生成 Voronoi 颗粒 的 奇异 性 。 

在 步 又 二 中 ,如 图 3 所 示 , 首 先 在 步骤 一 的 基础 上 
基于 所 有 随机 圆 的 圆心 生成 三 角形 网 格 , 再 以 三 角形 
的 名 接 圆 圆 心 作为 Voronoi 颗粒 的 角 点 生成 Voronoi 
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颗粒 群 。 由 此 得 到 的 Voronoi 颗粒 在 很 大 程度 上 继 
承 了 步骤 一 中 随机 圆 的 级 配 信 息 。 


&K|2. 子 域 随机 圆 的 生成 及 子 域 的 复制 


Fig.2 Random circle generation in sub-domain 


and its replication 


eo. 
ee 
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二 如 图 4(a) 所 示 , 在 步骤 三 中 ,连接 前 面 Al 、A2、 

I A4 子 域 的 中 心 ,截取 出 中 心 A0 区 域 , 并 对 AO 
进行 复制 拼接 , 直到 整个 区 域 的 大 小 满足 实际 

建 异 的 需要 。 截 取得 到 的 AO 区 域 模型 以 及 拼接 后 
的 重型 都 具有 边界 处 结构 的 周期 性 , 可 以 为 开展 代 
Ae VEU TH É 7C ( representative volume element, RVE ) 
提供 条 件 。 
asp 


(b) 周期 性 区 域 的 拼接 
到 4 周期 性 区 域 的 截取 与 拼接 
Fig.4 The cutting-out and assembling of periodic regions 
如 图 4(b) 所 示 , Voronoi 颗粒 的 区 域 大 小 满足 
了 建立 红色 圆圈 所 示 的 岩石 圆 盘 建 模 的 需要 。 为 避 
免 过 短 边 对 后 续 DDA 计算 的 影响 ,拼接 完成 后 还 对 
具有 过 短 边 的 Voronoi 颗粒 进行 了 局 部 重组 优化 。 


(a) 周期 性 区 域 的 截取 


图 3 基于 随机 圆 圆心 生成 Voronoi 颗粒 


Fig.3 Generating Voronoi particles based on random circle centers 


在 步骤 四 中 ,对 前 面 的 Voronoi 颗粒 拼接 模型 进 
行 边 界 或 孔洞 切割 ,生成 实际 建 模 所 需 的 外 轮廓 或 
内 部 孔洞 。 如 图 5 所 示 为 所 建立 的 带 有 一 定 宽 度 中 
心 预制 裂纹 的 岩石 圆 盘 模 型 ,其 中 的 预制 裂纹 由 程 
序 通过 孔洞 切割 算法 自动 生成 。 在 切 制 过 程 中 , 切 
制 边界 处 会 华 有 极 小 或 畸形 的 Voronoi 单元 生成 ,为 
避免 其 对 DDA 计算 精度 和 速度 产生 影响 ,还 将 在 不 
影响 边界 形态 的 基础 上 对 这 些 单元 进行 局 部 重组 优 
化 。 以 上 得 到 的 Voronoi 颗粒 单元 模型 ,最 后 在 步骤 
五 中 被 转化 生成 DDA 的 计算 模型 。 


5 ” 带 预制 裂纹 的 圆 盘 模 型 及 其 边界 优化 


Fig.5 Disc with a pre-existing crack and boundary optimization 
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2 圆 盘 径 向 压缩 破坏 模拟 


DDA 的 每 个 块 体 处 于 常 应 力 状态 。 无 论 采 
用 三 角形 ”| 或 Voronoi 单元 划分 , 子 块 体 单元 DDA 
方法 都 是 通过 引入 虚拟 节理 对 连续 岩石 模型 进行 离 
散 , 子 块 体 单元 间 的 虚拟 节理 赋 以 连续 岩石 的 强度 ， 
从 而 实现 连续 岩石 变形 的 模拟 。 当 子 块 体 单 元 的 应 
力 水 平 达 到 岩石 强度 时 ,临近 虚拟 节理 转变 为 真实 
节理 , 即 发 生 和 宕 石 的 开裂 破坏 。 其 中 岩石 的 拉 伸 和 
剪 切 开 裂 破坏 分 别 采 用 最 大 拉 应 力 和 莫 尔 -库仑 强 
度 准 则 。 


2.1 完整 圆 盘 


= 
三 图 6 为 完整 圆 盘 的 几何 模型 和 加 载 及 约束 条 
件 忆 圆 盘 通过 下 端 回 定 板 和 上 端 位 移 加 载 板 实 现 径 

Ati ,并 在 上 端 加 载 板 与 圆 盘 间 设 置 一 个 测量 板 。 


ji V=2 mm/s 


202307.00 
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R6 SER TANA 
Fig.6 An intact disc model 

测量 板 作为 一 个 常 应 力 DDA 抉 体 ,其 了 方向 应 

力 能 够 实时 反映 圆 盘 在 加 载 方 向 上 的 整体 应 力 水 

平 。 圆 盘 直 径 100 mm, 圆 盘 上 下 端 设置 中 心 角 为 

10° 的 平台 , 且 与 固定 板 及 测量 板 间 的 接触 面 设置 为 


. .se 
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光滑 面 ,从 而 缓解 加 载 端的 应 力 集中 。 加 载 板 向 下 
的 加 载 速度 取 为 2 mm/s, 圆 盘 单 元 数 考虑 在 约 1 500 ~ 
7 500 范 围 内 ,每 间隔 约 1 500 单元 数 建立 1 个 模型 ， 
从 而 对 圆 盘 的 辟 裂 破坏 模拟 进行 网 格 尺寸 敏感 性 分 
析 。 本 研究 暂 不 考虑 Voronoi 颗粒 单元 级 配对 圆 盘 
破裂 的 影响 ,因此 所 有 模型 都 采用 尺寸 相对 均匀 的 
网 格 。 圆 盘 的 力学 参数 和 DDA 计算 控制 参数 取 值 
如 表 1 所 示 。 


(a) 子 块 体 单元 数 1 475 (b) 子 块 体 单元 数 3 060 


(c) 子 块 体 单元 数 4 525 (d) 子 块 体 单元 数 6 001 


(f) 完整 岩石 圆 盘 实 验 结果 ™ 


(e) 子 块 体 单元 数 7 521 
图 7 不 同 子 块 体 单 元 数 下 的 完整 圆 盘 臂 裂 
的 DDA 模拟 结果 及 实验 结 


Fig.7 DDA simulation results with different number of 


sub-blocks and experimental result of an intact disc 


表 1 圆 盘 力 学 参数 与 DDA 计算 控制 参数 


Tab.1 Disc mechanical parameters and DDA computing control parameters 


弹性 模 量 /CPa ” 泊 松 比 。 抗 拉 强 度 AMPa ”内 摩擦 角 /(°) 


内 聚 力 /MPa 


弹 得 刚度 /GPa 时 间 步 长 /s 最 大 位 移 比 


20 0.2 1.0 35 


图 7 为 不 同 子 块 体 单元 数 完整 圆 盘 劈 裂 的 DDA 
模拟 结果 及 其 与 实验 结果 的 比较 。 可 以 看 出 ,不 同 
单元 数 条 件 下 圆 盘 的 最 终 臂 裂 破坏 形态 基本 一 致 。 
增加 子 块 体 单 元 数 可 以 使 初始 裂纹 更 接近 圆 盘 中 


5x10-* 


心 ,裂纹 贯穿 路 径 更 接近 理想 路 径 , 裂 纹 带 变 罕 LEY 
裂 路 径 更 清晰 ,从 而 与 实验 破坏 形态 更 加 接近 。 

图 8 是 子 块 体 单元 数 为 7 521 时 圆 盘 的 臂 裂 破 
坏 过 程 。 可 以 看 出 ,在 第 13 790 时 步 , 圆 盘 初 始 裂纹 


8.0 100 1x107 
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出 现在 接近 圆 盘 中 心 位 置 ,说 明 此 时 该 位 置 的 最 大 
主 应 力 率先 达到 材料 的 抗 拉 强 度 1.0 MPa。 随 着 加 
载 板 压 缩 位 移 的 继续 增加 ,裂纹 沿 加 载 方向 扩展 并 
贯穿 整个 圆 盘 ,从 而 发 生 臂 裂 破坏 。 需 要 说 明 的 是 ， 
为 了 较 好 地 捕捉 准 静 态 渐进 开裂 破坏 中 应 力 逐 渐 释 
放 对 开裂 破坏 过 程 的 影响 ” ,计算 中 1 个 时 步 内 仅 
aii din jmd sensi, a e dd 


e (c) 第 13 900 时 步 (d) 第 14 050 时 步 
图 8 子 块 体 单元 数 为 7521 时 的 完整 圆 盘 劈 裂 破坏 过 程 
8 Splitting process of intact disc with 7 521 sub-blocks 
ORI 为 不 同 子 块 体 单元 数 条 件 下 测量 板 的 了 方 
向 应 力 时 间 曲 线 。 
1.8 
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图 9 不 同 子 块 体 单元 数 下 完整 圆 盘 的 了 方向 整体 应 力 


Fig.9 The integral stress-time curve in Y direction of the 


0.012 0.015 


intact disc with different number of sub-blocks 


可 以 看 出 , 圆 盘 的 了 方向 整体 应 力 先 呈 线性 增 


加 。 随 着 加 载 的 持续 , 圆 盘 中 心 附近 的 子 块 体 单 元 
应 力 首先 达到 岩石 的 拉 伸 强度 ,从 而 出 现 初始 裂纹 
并 沿 加 载 方向 扩展 ,并 在 较 小 时 间 段 内 完成 圆 盘 的 
臂 裂 破坏 , 圆 盘 的 了 方向 整体 应 力 随 之 迅速 下 降 。 
子 块 体 单元 数 越 小 时 , 圆 盘 表现 出 更 高 的 整体 强度 ， 
这 是 因为 一 方面 单元 数 的 减少 将 降低 了 应 力 场 描述 
的 准确 性 , 另 一 方面 更 少数 量 的 虚拟 裂纹 为 裂纹 的 
发 育 提供 了 更 少 的 路 径 机 会 

根据 圆 盘 内 初始 裂纹 出 现时 测量 板 的 了 方向 应 
力 可 以 换算 得 到 此 时 圆 盘 的 径 向 外 载荷 已 ,并 进而 

根据 公式 ox = 2P/«D (D 为 圆 盘 直径 ) 计算 得 到 此 

时 圆 盘 中 心 的 民 方 向 应 力 。 表 2 列 出 了 不 同 子 块 体 
单元 数 条 件 下 圆 盘 内 初始 裂纹 出 现 的 时 间 以 及 此 时 
的 径 向 加 载 力 及 所 求 圆 盘 中 心 应 力 。 可 以 看 出 : 子 
块 体 数 增多 时 , 出现 初 始 裂纹 时 圆 盘 中 心 式 方向 应 
力 总 体 上 越 趋 近 于 1.0 MPa ,计算 的 精度 更 高 ; 当 子 
块 体 数 为 6 001 及 以 上 时 ,计算 误差 均 在 1. 596 以 
内 ,达到 了 可 接受 的 计算 精度 。 

表 2 不同 子 块 体 数 下 圆 盘 开裂 时 间 点 与 对 应 载荷 


Tab.2 The crack initiation time of disc with different number 


of sub-blocks and the corresponding load 


子 块 体 单元 数 1475 3 060 4525 6001 7521 
起 裂 时 间 /s | 0.01323 0.01340 0.01367 0.01357 0.01378 
外 载荷 /kN 166.2 161.0 160.6 157.9 159.3 

中 心 wx /MPa 1.058 1.025 1.022 1.005 1.014 


2.2 45?" 倾 角 预 制 裂纹 圆 盘 


图 10 为 一 个 带 a =45° 倾 角 中 心 预制 裂纹 的 圆 
盘 模型 。 圆 盘 直 径 为 100 mm ,预制 裂纹 的 长 度 与 宽 
度 分 别 为 30 mm 和 0.25 mm ,与 文献 [2] 相 应 实验 中 
的 设置 一 致 ,加 载 板 向 下 的 加 载 速度 取 为 0.5 mms, 
圆 盘 的 力学 参数 和 DDA 计算 控制 参数 与 2. 1 节 一 
致 ,并 仍然 在 大 约 1500 ~7 500 范围 内 每 隔 约 1 500 
单元 数 开 展 模拟 分 析 。 
图 11 为 不 同 子 块 体 单 元 数 下 圆 盘 破坏 的 模拟 
结果 及 其 与 实验 结果 的 比较 。 不 同 子 块 体 单元 数 条 
件 下 圆 盘 的 起 裂 点 都 是 在 预制 裂纹 的 两 个 尖端 处 ， 
之 后 裂纹 扩展 方向 逐步 趋 近 于 加 载 方向 并 最 终 贯穿 
整个 圆 盘 。 子 块 体 单元 数 越 多 ,对 应 的 扩展 路 径 越 
清晰 ,裂纹 带宽 度 越 罕 ,破坏 路 径 越 接近 实验 结果 。 
从 图 12 中 的 了 方向 整体 应 力 时 间 曲 线 可 以 看 出 , 当 
子 块 体 单元 数 较 少 时 , 圆 盘 整体 强度 相对 更 高 。 当 
子 块 体 数 增加 到 约 4 500 及 以 上 时 ,曲线 基本 收 和 敛 , 说 
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明 此 时 单元 数 对 计算 结果 的 影响 已 经 很 小 。 此 外 , 当 


单元 数 从 3 082 增加 到 4 553 时 ,曲线 出 现 了 跳 变 ,这 


应 该 是 本 组 算 例 中 所 取 的 单元 数 间隔 偏 大 所 致 。 


图 10 带 45" 倾 角 预 制 裂纹 的 圆 盘 模型 
Fig. 10 Disc model with a 45? pre-existing crack 


(b) 子 块 体 单元 数 3 082 


(a) 子 块 体 单元 数 1 492 


(d) 子 块 体 单元 数 6 027 


(c) 子 块 体 单元 数 4 553 


(e) 子 块 体 单元 数 7 540 (f) 实验 结果 ” 
图 11 不 同 子 块 体 单元 数 下 45" 预 制 裂纹 
圆 盘 破坏 模拟 结果 和 实验 结果 
Fig.11 DDA simulation results with different number of 


sub-blocks and experimental result of the disc 


with a 45? pre-existing crack 


D 
周 大 波 ， 等 ;岩石 圆 盘 径 向 压缩 破坏 的 Voronoi 子 块 体 单元 ERORA (FERT 29 
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图 12 不 同 子 块 体 单元 数 下 45" 预 制 裂 纹 
圆 盘 的 了 方向 整体 应 力 


Fig. 12 The integral tress-time curve in Y direction of the disc 


with 45? pre-existing crack under different number of sub-blocks 


2.3 不 同 倾角 预制 裂纹 圆 盘 


在 2.2 节 的 基础 上 ,本 节 考 虑 圆 盘 中 心 预制 裂 
纹 倾 角 a 由 0° 到 75° 变 化 的 情况 ,每 隔 15° 进 行 模拟 
分 析 , 圆 盘 的 子 块 体 单元 数 均 在 7 500 左右 。 图 13 
和 图 14 分 别 为 不 同 倾角 中 心 预制 裂纹 圆 盘 破坏 的 
模拟 结果 和 实验 结 


(a) a=0° (b) a=15° 


(d) a=45° 


(e) a=60° 
13 不同 倾 角 中 心 预制 裂纹 圆 盘 的 模拟 结 


Fig.13 Simulation results of discs with different 


(f) a=75° 


pre-existing crack inclinations 
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图 14 不 同 倾角 中 心 预 制 裂纹 圆 盘 的 实验 结果 
Fig.14 Experimental results of discs with different 
[2] 


pre-existing crack inclinations 

一 可 以 看 出 , 除 a =75° 的 情况 外 ,DDA 模拟 得 到 
的 图 盘 开裂 破坏 路 径 与 相应 的 实验 结果 高 度 吻合 。 
HENE 0° -60* 范 围 内 时 , 圆 盘 的 起 裂 点 都 是 在 预制 
发 妆 的 两 个 尖端 处 , 随后 裂纹 扩展 路 径 逐 渐 趋 近 于 
垢 惰 方 向 并 最 终 贯穿 整个 圆 盘 。 当 a =75° 时 , DDA 
模 私 和 实验 中 的 初始 裂纹 都 不 是 出 现在 预制 裂纹 尖 
HEDDA 模拟 中 初始 裂纹 出 现在 离 预制 裂纹 两 个 
类 油 一 定 距离 位 置 的 初始 裂纹 面 处 ,这 与 文献 [22] 

实验 现象 是 一 致 的 。 而 在 图 14 的 实验 中 ,初始 
要 融 出 现在 圆 盘 中 心 位 置 ,这 是 因为 受到 了 圆 盘 中 
心 独 始 加 工 孔 的 影响 。 

图 15 是 a 为 60 和 75。 时 ,考虑 圆 盘 中 心 直 径 
Hmm 的 小 孔 时 的 模拟 结果 。 可 以 看 出 DDA 模拟 
结 时 与 实验 结果 吻合 ,证 明 在 a = 60* 时 , 圆 孔 对 开 
裂 尾 坏 影响 可 以 忽略 。 事 实 上 ,更 多 的 计算 表明 , 当 
a <60° 时 , 圆 筷 对 圆 盘 开裂 破坏 路 径 的 影响 都 可 以 
忽略 。 而 a =75° 时 , 圆 孔 附近 的 应 力 集中 导致 了 初 
始 裂纹 从 圆 盘 中 心 发 育 。 以 上 结果 表明 DDA 模拟 
很 好 地 捕捉 到 了 圆 盘 受 载 过 程 中 的 变形 和 开裂 破坏 
特征 。 


多 15 带 中心 小 孔 预 制 裂 纹 圆 盘 的 模拟 结果 


Fig. 15 Simulation results of pre-existing crack 


discs with a center hole 


ChinaXiv& ERBTI a 4o % 
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本 研究 提出 了 一 种 基于 随机 圆 的 Voronoi 单元 
模型 建立 方法 ,该 方法 能 够 考虑 颗粒 的 级 配 特性 , 具 
有 较 高 的 生成 效率 ,并 能 保证 必须 要 的 周期 性 。 径 
向 压缩 加 载 下 完整 圆 盘 和 和 带 预 制 裂纹 圆 盘 开裂 破 
坏 的 子 块 体 单元 DDA 方法 结果 模拟 表明 : 较 小 的 子 
块 体 单元 数 将 得 到 更 高 的 圆 盘整 体 强度 ; 随 着 子 块 
体 单 元 数 的 增 大 , 圆 盘 的 整体 响应 曲线 趋 于 收敛, 破 
坏 起 裂 位 置 更 加 接近 真实 ,开裂 破坏 路 径 更 加 清晰 ; 
采用 较 大 的 单元 数 将 能 够 得 到 与 实验 结果 一 致 的 开 
裂 破 坏 模拟 结果 ,并 且 模 拟 还 准确 的 再 现 了 圆 盘 中 
心 小 孔 对 开裂 破坏 路 径 的 影响 。 

使 用 Voronoi 子 块 体 单元 DDA 方法 能 够 有 效 模 
拟 岩 石 的 开裂 破 坏 过 程 。 下 一 步 将 开展 Voronoi W 
粒 级 配 及 形态 对 岩石 开裂 破坏 模拟 的 影响 研究 。 
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